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The sizes of the field galaxies with / < 25 have been measured in the NTT Deep Field. 

^ ■ Intrinsic sizes have been obtained after deconvolution of the PSF with a multigaussian method. 

The reliability of the method has been tested using both simulated data and HST observations 

0^ ' of the same field. The distribution of the half light radii is peaked at rhi — 0.3 arcsec, in good 

^^ ■ agreement with that derived from HST images at the same magnitude. An approximate morpho- 

logical classification has been obtained using the asymmetry and concentration parameters. The 
^ ' intrinsic sizes of the galaxies are shown as a function of their redshifts and absolute magnitudes 

r^ . using photometric redshifts derived from the multicolor catalog. While the brighter galaxies 

with morphological parameters typical of the normal spirals show a flat distribution in the range 

f^ ■ r^ = 1 — 6 kpc, the fainter population at 0.4 < z < 0.8 dominates at small sizes. To explore the 

significance of this behaviour, an analytical rendition of the standard CDM model for the disc size 
evolution has been computed. The model showing the best fit to the local luminosity function 
and the TuUy-Fisher relation is able to reproduce at intermediate redshifts a size distribution 
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f^ , in general agreement with the observations, although it tends to underestimate the number of 

I ' galaxies fainter than Mb '^ — 19 with disk sizes r^; ^ 1 — 2 kpc. 

C^ ■ Subject headings: galaxies: fundamental parameters - galaxies: evolution - galaxies: formation 
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1. INTRODUCTION 

Deep multicolor surveys from HST and ground based telescopes coupled with deep spectroscopy on ten 
meter class telescopes (Keck, VLT) are leading to a remarkable progress in our understanding of the galaxy 
formation and evolution. In particular, there has been great improvement in the knowledge of the epoch of 
star formation. The first determinations of the global star formation history suggested a peak at z « 1.5 
with a further decline from z ~ 2 to z « 4 (Madau et al. 1996; Connolly et al. 1997). In this scenario, more 
than 80% of the stars would have been formed recently at z < 2. 

Such a picture would be in rough agreement with the cold dark matter (CDM) models of galaxy for- 
mation, which predict the star formation to be a gradual process governed by gas cooling and Supernovae 
feedback. The timescale of the process is strongly related to the rate at which the mass is hierarchically 
assembled to form larger galaxies from smaller structures (White & Frenk 1991; Cole et al. 1994; Kauffmann, 
White & Guiderdoni 1993; Baugh et al. 1998). 

However, the above scenario is presently challenged by recent indications that the observed decline at z > 
2 is affected by incompleteness in the color selection and by dust extinction. In fact, more recent observational 



evidence limited to the optically selected galaxies, points towards a flat SFR for z > 2 (Giallongo et al. 1998, 
Pascarelle et al. 1998, Steidel et al. 1999). 

A more precise test of the galaxy formation scenarios involves the study of the size and morphological 
distribution of galaxies with cosmic time, since this is related to the main feature of the hierarchical models, 
namely, the growth of mass and sizes of typical condensations. Indeed, although the global star formation 
history is a crucial issue in galaxy formation, it is an integrated quantity; even a flat SFR could be consistent 
with the CDM picture by tuning the parameters describing the gas cooling and the Supernovae feedback in 
galaxies. 

The first observational investigations of the differential evolution of the galaxy population over the 
redshift range < z < 1 were carried out by Lilly et al. (1995), Cowie et al. (1996) and Ellis et al. 
(1996) by means of spectroscopic redshift surveys. The observed decrease with the cosmic time of the blue 
star-forming galaxy population showing strong emission lines has been further investigated using the Hubble 
Space Telescope. The decline is a mixture of stellar evolution in normal spirals and strong evolution in low 
luminosity irregulars (Schade et al. 1995, Brinchmann et al. 1998, Lilly et al. 1998) whose local counterparts 
now comprise the faint end of the local luminosity function (Zucca et al. 1997). 

Subsequently, the study of the Hubble Deep Field (Williams et al. 1996) has resulted in a strong 
improvement in the knowledge of the galaxy morphology at high redshift (Abraham et al. 1996, Odewahn et 
al. 1996, Driver et al. 1998) showing a high fraction of faint irregular galaxies. Most of these galaxies have 
small sizes resulting in a physical extent of the order of 3 — 5 kpc at z ^ 1 (Glazebrook et al. 1998). Such 
small sizes, if compared with the present day galaxy morphology included in the Hubble sequence, could 
imply physical growth and merging occurring with increasing cosmic time (Bouwens, Broadhurst and Silk 
1997), in qualitative agreement with the hierarchical clustering picture. On the other hand, the indications 
for a slow evolution of the spiral population constituted by the presence of large, relatively bright spirals at 
z ^ 0.6 (Schade et al. 1995) are difficult to fit in the picture of the CDM-like models. Indeed, the dissipative 
formation of galactic disks in CDM halos predicted by the hierarchical clustering scenario produces spirals 
whose size strongly decreases as the formation redshift increases. For this reason, gas cooling in the CDM 
halos has to be delayed to low redshifts (z < 1). The comparison becomes more stringent for the lower 
luminosity spirals where the predicted sizes are very small at high redshifts. 

A wide field survey of faint galaxies observed with the Hubble Space Telescope would be the ideal 
sample to clarify the issue and to provide a differential and quantitative test to CDM galaxy formation 
theories. However deep photometry from HST can only be obtained on small areas. On the contrary, wide 
field surveys of faint galaxies can be obtained using ground-based instrumentation but at lower resolution 
(FWHM- 0.5 - 1 arcsec). 

In this paper we have exploited the deep photometric information of the NTT Deep Field (Arnouts et 
al. 1999) obtained under condition of subarcsec seeing at the ESO NTT telescope to derive morphological 
information after appropriate seeing deconvolution (Sect. 2). The reliability of the deconvolution technique 
has been tested by means of numerical simulations and particularly by checking on somewhat shallower HST 
images of the same galaxy field retrived from the HST archive. The results show that it is possible to derive 
the radial intensity profiles of galaxies with / < 25 from ground-based imaging with FWHM< 1 arcsec (Sect. 
3). 

The derived intrinsic angular sizes are then converted into physical sizes adopting reliable photometric 
redshifts estimated from the UBVRIJK multicolor distribution of each galaxy in the catalog. The galaxy 
morphological distribution in redshift and luminosity is shown in Sect. 4 and is compared with that predicted 



by hicrachical clustering models for the disk formation in Sect. 5. The conclusions are presented in Sect. 6. 



2. THE GALAXY CATALOG 

The galaxy catalog has been obtained from the NTT Deep Field (NDF) (Arnouts et al. 1999). The 
optical observations were obtained in the Johnson BV, and Gunn R,I broadband filters with the Ik x Ifc 
SUSI camera on the ESO New Technology Telescope. The resulting field of view is ~ 5 arcmin^. 

The 3(7 magnitude limits arc B ~ 27.5, V = 27, R = 26.5, / = 26.5, the average seeing was ~ 0.8 arcsec 
(see Arnouts et al. 1999 for further details). 

Deep ultraviolet images have been added with the SUSI-2 camera at the NTT reaching a 3(t upper limit 
oi U = 26. Also deep near infrared observations have been added with the SOFT instrument at the NTT 
reaching 3a magnitude limits of J ~ 23.5 and K' ~ 21.7. 

The galaxy catalog has been extracted from the sum of the BVRI weighted frames to detect objects at 
very faint level in the four bands (Arnouts et al. 1999). The object detection and photometry have been 
performed with the Sextractor image analysis package (Berlin & Arnouts 1996). Total magnitudes in the 
UBVRIJK bands were computed from isophotal magnitudes for galaxies brighter than 1=23.25 and from 
aperture magnitudes in 2.2 arcsec corrected to 5 arcsec for fainter objects. Colors are computed in a 2 arcsec 
diameter aperture in the same sky area in all the bands after correction for the different seeing conditions. 
Details on the data reduction and photometry are fully discussed in Fontana et al. (1999). 

To each galaxy in the UBVRIJK catalog, a photometric redshift estimate has been assigned with the 
same best fitting procedure applied in the contiguous BR1202-0725 field (Giallongo et al. 1998). This 
has been obtained through a comparison of the observed colors with those predicted by spectral synthesis 
(GISSEL) models (Bruzual & Chariot 1993) including UV absorption by the intergalactic medium and dust 
reddening. The catalog of the photometric redshifts is presented in Fontana et al. (1999). The resulting 
typical redshift accuracy is Az ^ 0.1 up to z ~ 1.5 and Az ~ 0.15 at larger redshifts. 

The redshift distribution in the NTT Deep Field is shown in Fig. 1 for galaxies with / < 25. As derived 
in other fields of similar depth, the bulk of the galaxies present are at intermediate redshifts z ~ 0.5 — 0.8 
with a tail in the distribution up to z ^ 4. 

The shape of the optical counts in this field suggests negligible correction for incompleteness at / = 25 
(Arnouts et al. 1999) thus the morphological study described in the next sections should be considered 
complete to this magnitude limit. 



3. THE RADIAL INTENSITY PROFILES OF FAINT GALAXIES 

The estimate of the intrinsic radial intensity profiles derived from deep ground based images of faint 
galaxies is a complex task. Indeed, HST images of galaxies with / ~^ 25 have shown that most of them have 
typical sizes corresponding to half light radii r^i ^ 0.3 arcsec which are smaller than the typical HWHM^ 0.5 
arcsec present in ground based deep images. 

Accurate and reliable deconvolution techniques must be applied to obtain any statistical information 
about the intrinsic half light radii of faint galaxies observed from the ground. Any deconvolution technique 



depends of course on the accurate knowledge of the instrumental profile (PSF) and on the intrinsic shape 
assumed for the galaxies. An accurate reconstruction of the instrumental profile is generally done through 
an average of the stellar profiles with high s/n ratio in images with good oversampling. 

Concerning the intrinsic shape of the galaxies, instead of assuming a specified profile, we adopted the 
multigaussian deconvolution analysis studied in detail by Bendinelli (f 99f ) and also recently applied to faint 
elliptical galaxies in the HDF by Fasano et al. (1998). 

For each object detected in the sum of the BVRI weighted frames, a radial intensity profile was obtained, 
sampling the light into elliptical annuli with axes proportional to the object intensity-weighted second order 
moments. Following Fasano et al. (1998) we chose to carry out the analysis on a circularly symmetrized 
profile, where the radial distance along the semi-major axis a„ is r„ = %, where e is the Sextractor elongation 
parameter. 

A multi gaussian expansion of these symmetrized profiles has been developed as described below. If I{r) 

is the true, deconvolved profile of the object, it is possible to obtain a simple but effective approximation of 

it: 
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where L^ is the total luminosity of the source, Sj and bj the spread parameter and the weight of each 
gaussian, respectively. We adopt M = 6 as a compromise between the need of sampling the profile with the 
maximum allowed number of gaussians and the need of performing a well-posed fit. Weighting parameters 
are obviously not independent: their sum must be equal to one. 

In an analogous way an expansion of the normalized point spread function in the frame was carried out, 
so that it was possible to have an analytical solution of the convolution integral: 
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where the terms with index "i" refer to the multi gaussian expansion of the PSF. Once the PSF expansion 
is known, it is possible to fit Eq. to the observed points, obtaining the spread and weighting parameters 
Sj, bj of the deconvolved profile. 

In order to verify the reliability of the deconvolution technique in our particular case, we have performed 
two checks. First of all we have produced a set of simulations specifically designed to reproduce typical 
conditions in our sample. Second, we have compared, for a bright galaxy subsample, intrinsic half light radii 
derived from the I band galaxy profiles taken at the NTT with that derived from a relatively deep image in 
the I band taken at the WFPC2 on HST. 

In the first check, intensity profiles have been reproduced as in the observed NTT deep field assuming the 
same pixel sampling and the same average seeing (FWHM^ 0.8 arcsec). Assuming an intrinsic exponential 
profile, a series of synthetic images were constructed using different half light radii ranging from r^i = 0.1" 
to rhi = 0.9" with a step of 0.1". An average of 25 random objects were computed for each radius, assuming 
a total magnitude of / ~ 23.5 and of J ~ 24.7 which are typical values for a bright and a faint object in 
our sample. In this way the s/n ratio is constant in the simulations. In order to reproduce the average 
seeing conditions in the NDF, the background subtracted image of a bright star was selected in the field. 
Its normalized profile was then convolved with the synthetic images of the disk galaxies. The convolved 



two dimensional profiles were randomly inserted in regions of the NTT deep field far from very bright 
objects to reproduce with the appropriate pixel size and noise levels the observed NTT galaxies. Finally, the 
multigaussian deconvolution technique was applied to the synthetic data. We first notice that we have no 
selection bias against galaxies with large size (r/j; ^ 1 arcsec) and low surface brightness down to / ~ 24.7 
since all the synthetic objects were detected. 

The results are shown in Fig. 2, where the error bars represent the dispersion around the mean due 
to noise in the background subtraction. A good match between the intrinsic and measured half light radii 
has been obtained up to rhi ~ 0.7 arcsec. For larger values, a slight underestimate in the measured values 
appears particularly for the faint subsample. In any case it can be seen that, even for the faint galaxies with 
I ^ 25, the overall correlation between intrinsic and measured half-light radii is preserved in such a way that 
an intrinsically large, faint object, e.g. with rhi ~ 0.7", can not be detected as a small sized one, e.g. with 
Vfii ^ 0.1". The simulation shows that the fraction of small size galaxies observed in the NDF is real and is 
not due to intrinsically larger objects which have been shrunk by noise effects. 

For the brighter galaxies we had the opportunity to check directly the reliability of our deconvolution 
technique comparing our I band image photometry of the NTT deep field with the I band WFPC2 image 
photometry taken from the HST archive. 

We have applied to the WFPC2 galaxies the same deconvolution technique using an average PSF derived 
from the available stellar profiles in the field. Since the WFPC2 image is not as deep as the NTT I band 
image and the overlap between the two images is limited to less than half of the field, the comparison 
was limited to a few relatively bright objects with / < 24.5. In Fig. 3 the deconvolved half light radii 
derived from ground based and HST I band images are compared. The agreement is very good with a 
rms difference of 0.05. A similar agreement for the galaxy scale lengths has been found by Schade et al. 
(1996) for a sample of CFRS galaxies two mag brighter, Iab < 22.5. Thus, although HST images allow 
more detailed morphological studies with a decomposition of the intensity profiles into more components 
(e.g. bulge plus disk decomposition), the statistical information is limited by the small field of view. Large 
statistical information on the galaxy rhi from wide field ground based images can be used to constrain the 
models of galaxy formation. As a first step, we have applied this technique to the multicolor NTT deep field. 
The results are discussed in the next Section. 



4. THE DISTRIBUTION OF THE GALAXY SIZES IN LUMINOSITY AND REDSHIFT 

The distribution of the half light radii rhi for all the galaxies in the multicolor catalog with / < 25 
is shown in Fig. 4. The median value of the distribution is rhi — 0.3", which confirms previous findings 
from the HST medium deep survey (MDS). In fact, the derived distribution is very similar to that of the 
MDS (Roche et al. 1998). This is an important result since it has been obtained by means of ground based 
observations. 

Assuming that all the galaxies are characterized by a disk with a given scale length r^ — r/ij/1.68, we 
have computed the linear size for each galaxy as a function of redshift using the color estimated redshifts 
as discussed in section 2. The resulting luminosity-size relation is shown in Fig. 5 in two redshift intervals. 
Uncertainties in the redshift estimates Az ^ 0.1 produce uncertainties in the sizes and luminosities ^ 5% 
and ~ 30%, respectively. 

Fig. 5 shows a concentration of faint galaxies with r^ < 3 kpc that is present at Mb > —19. This is 



more evident in the 0.4 < z < 0.8 interval where the galaxy distribution extends down to Mb ^ —16. 

The histogram of the size distribution is shown in Fig. 6 in the same redshift intervals where a peak of 
low-luminosity dwarf galaxies appears for r^; < 2 kpc. 

However the assumption of a disk-like structure for the majority of the faint galaxies can no longer be 
valid because of the presence of a noticeable fraction of dwarf galaxies with irregular morphology. Although 
it is not clear if the images of the high z galaxies seen at progressively shorter wavelengths can produce 
a spuriously increasing number of irregular galaxies, we have attemped to exclude galaxies with irregular 
morphology following the method described in Abraham et al. (1996). With this method, galaxies are 
separated into three classes (irregulars, spirals, SO-ellipticals) according to the values of their asymmetry 
and concentration parameters. The asymmetry parameter is defined as: 

1 VIZ- ■ -/i^"l 

where / and J^^" are the background-subtracted intensity of the original image and the same after 180-degree 
rotation around the baricenter, respectively. Here the sum is extended to the isophotal detection area of 
every source, and k takes account of the bias introduced in A by the underlying portion of the sky with the 
same area and shape as the object. For the concentration parameter we have: 

C„ = ^Mf£^ (4) 

where S refers to a circular aperture with a radius r so that it has the same number of pixels as the isophotal 
detection area of the object. In an analogous way Sa is an inner circular area with radius r^ — ar. We 
chose to adopt a = 0.4 instead of a = 0.3 because we found that this larger value takes account of the 
seeing blurring of the sources in the NDF compared to the HST high resolution images. The distribution in 
the asymmetry-concentration plane shown in Fig. 7 is in good agreement with the Abraham et al. (1996) 
distribution derived from the Hubble Deep Field galaxies of comparable magnitude. 

In Fig. 7 regular spiral galaxies are distributed in the central region delimited by the two straight lines 
defined in Abraham et al. (1996). Their distribution in the luminosity-size plane is also shown by filled 
symbols in Fig. 5. We have reproduced in Fig. 8 the size distribution of the normal spirals in the same 
redshift interval as in Fig. 6. It can be seen that even reducing the sample to the more regular spirals, the 
size distribution in Fig. 8 is very similar in shape to that of the whole sample shown in Fig. 6 with a peak 
at very small sizes r^ ~ 1.5 kpc. This is expected since the distribution of disk-dominated systems in the 
luminosity-size plane (filled circles in Fig. 5) follows the distribution of the whole sample. 

This behaviour has interesting consequences for cosmological scenarios of galaxies formation. 



5. EVOLUTION OF GALACTIC DISKS IN CDM COSMOLOGY 

The size distribution of different spectral types as a function of redshift is in principle a powerful tool 
to constrain models of galaxy formation in hierarchical CDM cosmogonies. 

Indeed, in these models the galaxy formation proceeds through subsequent merging of small units into 
larger ones, corresponding to a hierarchical growth of the typical mass of virialized condensations. Since 
these models predict galaxies in the range 1 < z < 3 to be still assembling, the corresponding size growth 
should be detectable by present-day surveys. 



Quantitatively, to compare the predictions of the CDM models with deep imaging observations, one 
needs to express the outcomes of the hierarchical clustering picture in terms of observable properties, i.e., 
luminosities and sizes. The first goal is usually pursued through the use of semi-analytical models (White & 
Frenk 1991; Cole et al. 1994; Kauffmann, White & Guiderdoni 1993; Baugh et al. 1998), which are able to 
predict luminosity functions and TuUy-Fisher relations in reasonable agreement with existing observations. 
As for the galaxy sizes, quantitative attempts to relate them to the mass growth of cosmic structures have 
been done only for the regular spiral population (Mo, Mao & White 1998). 

To compare with observations we have to connect the galaxy sizes with their luminous properties. To 
this aim, we first used the Mo et al. (1998) model to connect the scale length of any exponential disk Rd to 
the dark matter (DM hereafter) circular velocity Vc', then we used the approach of semi-analytical models 
to connect Vc to the luminous properties of galaxies. The predicted relations among luminosity, size and 
redshift of the field spirals have then been used to derive statistical distributions which have been compared 
with the observed ones. 



5.1. Size-circular velocity relation 

In the framework of present theories of galaxy formation the size of the disk R^ can be related to the 
circular velocity of the disk (vd) and of the DM {v^) halo, assuming that the mass and angular momentum 
of the baryonic material which settles into the disk are fixed fractions - rrid and j^ - of the DM mass m 
and angular momentum J (see Mo, Mao & White 1998). The latter can be expressed in dimensionless units 
A = J/(|i?|^/^/GTO^/^), whose distribution p(A) is known from several N-body experiments (see Warren et 
al. 1992; Cole & Lacey 1996; Steinmetz & Bartelmann 1995). 

For a galaxy DM density profile of the form given by Navarro, Frenk & White (1996, NFW hereafter), 
it is found that 
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where H{z) is the Hubble constant at redshift z and c = rs/r2oo is the concentration factor expressed in 
terms of the scale radius r^ in the NFW distribution and of the virial radius r2oo • The function fc is related 
to the total energy of the NFW profile while fn includes also the gravitational effect of the disk in the density 
distribution. Fitting formulas for the last two functions are provided by Mo et al. (1998). 

Note that in the above model the disk contributes to the rotation velocity, increasing by more than 20 
% the value of Vc given by the DM halo. We have taken this effect into account in our computations. The 
power spectrum of the initial DM density perturbation (from which the DM halos are thought to be formed 
via gravitational instability) enters the computation through the concentration factor c. 

According to the above discussion, and following the lines in Mo et al. (1998), we have computed the 
Rd — Vd and the Rd — Vc relations for different CDM perturbation spectra. In order to compare the results 
with the outcome of the data analysis we have converted the Rd — Vd (or Rd — Wc) relationships into Rd — Mb 
as described in the next section. 



5.2. Size-magnitude relation using a model for galaxy formation 

To implement the present understanding of the process of galaxy formation and evolution - and thus to 
associate magnitutes to the DM galaxy masses - we propose an analytical rendition of the semi-analytical 
approaches introduced by Kauffman, White & Guiderdoni (1993) and by Cole et al. (1994) and subsequently 
developed by several authors (see Somerville & Primack 1998). In these models the DM halos form through 
the gravitational instability of initial random density fluctuations whose power spectrum depends on the 
nature of the DM. The baryons within the DM halo, heated to the virial temperature of the halo, are assumed 
to be distributed in three phases: a cold phase where the gas is radiatively cooled and star formation may 
occur, a condensed phase (stars) and a hot phase (gas re-heated by Supernovae explosions). The hierarchical 
merging of the halos (followed by the possible merging of the galaxies contained in the halos) and the 
luminosity evolution due to the evolution of the stellar population are also included in the model. 

The interplay among these phases and the details of our implementation of the hierarchical models of 
galaxy formations are given in Appendix A. Following the lines in the literature, the parameters regulating 
the star formation have been chosen so as to match the local observed properties of the galaxies, i.e., the I- 
band TuUy-Fisher relation and the B-band luminosity function. However, since the disk velocity v^ is ^ 20% 
higher than that of the DM, matching the predictions of the model with both the above observables is critical. 
E.g., setting the star formation parameters to maximize the match to the observed TuUy-Fisher relation for 
a "Milky Way" galaxy with Vc~22Q km/s, yields an excess of bright galaxies in the local luminosity function 
(as noted by Somerville & Primack 1998). A further and complementary constraint is provided by the faint 
end of the luminosity function; minimizing the mismatch between the observed and predicted TuUy-Fisher 
relation at faint magnitudes leads to a luminosity function too steep at the faint end. 

For these reasons, we have selected the set of free parameters that minimizes the discrepancy of the 
model results with respect to both the local TuUy-Fisher relation and the local B-band luminosity function. 
For the angular momentum parameters (determining the Vc/vd ratio and thus affecting the predicted TuUy- 
Fisher relation) we have obtained jd = irid — 0.3. For the star formation parameters we have obtained the 
set reported in Appendix A. For the cosmology we have adopted a tilted CDM cosmogony (normalized to 
COBE with r^o = 1 and Hq ~ 50 km s^^ Mpc^^); the predictions of semi-analytical models with such a 
power spectrum are very close to the mixed (hot+cold) DM models, which are in fair agreement with the 
data for a wide set of observations (see Primack et al. 1995). With this recipe we have obtained a reasonable 
agreement of the model prediction with the local luminosity function and TuUy-Fisher relation (Fig. 9). The 
growth of the disk size with cosmic time, a characteristic prediction of CDM models, is shown in fig. 10. 
The size distributions, plotted at different z for two ranges of absolute magnitudes, are computed convolving 
eq. 5 over the full distribution p{X) of the halo angular momentum A that a galaxy with given Rd and Vd 
can have at any z. We have taken for p{\) the expression given in Mo et al. (1998). 



5.3. Results and comparison with the observations 

In Fig. 11 we show how the galaxy sizes are distributed in absolute magnitude in four redshift intervals. 
The dots are the observed data while the shaded regions correspond to the CDM predictions. The latter 
depend on the value of the dimensionless angular momentum A which is distributed according to the function 
p(A) in the interval 0.025 < A < 0.1, the extreme values being the 10% and 90% points of p(A), respectively. 
The solid line corresponds to A = 0.05, i.e., to the 50% point of the distribution p(A). 

Note that at bright {Mb < —19) magnitudes the observed points are approximately equally distributed 



above and below the solid line, and fall in the region allowed by the CDM model. In fact, for the bright 
galaxies (where a measurement of the disk rotation velocity has been possible) the model by Mo et al. (1998), 
relating size and circular velocity, is in agreement with the observations; since our model TuUy-Fisher relation 
relating v^ to Mg is also in agreement with observations (see Fig. 9), the matching between our predicted 
and observed distributions provides a consistency check of the theoretical approach used to derive the galaxy 
properties in CDM models. However, at fainter magnitudes (Mb > —19), the observed sizes tend to occupy 
systematically the small size region below the 50% locus of the angular momentum distribution. 

To investigate more quantitatively this point, the differential size distribution of galaxies with / < 25 
is shown in Fig. 8, for different redshift bins. Here the model results include the full p{\) distribution to 
take account of all the values of angular momentum that a galaxy with given R^ and Vd can have at a given 
z. To compare the model with the observed data in our sample, we have also introduced a magnitude limit 
at / < 25. This limit has the effect of selecting brighter (and hence larger) objects at higher z. This effect 
tends to balance the decrease with z of the disk size predicted by CDM models (as shown in fig. 10 at given 
absolute magnitude), resulting in an apparent lack of evolution with z of the distributions predicted by the 
CDM model (see Fig. 8). 

The comparison with the data shows a general agreement, although an excess of small-size objects is 
evident with respect to the CDM predictions. This points towards a galaxy size evolution from high to low 
z slower than predicted by CDM models, with the persistence of small size objects dominating in number at 
faint magnitudes and at intermediate to low redshifts. The differential distributions in Fig. 8 are plotted in 
terms of absolute numbers in our 5 arcmin^ field, so that the excess at small sizes over the CDM predictions 
is not affected by any possible incompleteness at the large/bright end of the distribution, which could be 
present at lower z in our limited sized field. 

We also note that the excess is robust with respect to any possible systematic errors in the selection of 
disk galaxies using the concentration-asymmetry plane shown in Fig. 7. Indeed, when the whole sample is 
considered (including also bulge and irregular systems) the resulting distributions shown in Fig. 6 are very 
similar to the distributions of the "spiral" sample shown in Fig. 8. In particular, the excess with respect to 
the CDM predictions at small sizes is still present. 



6. CONCLUSIONS 

The size distribution in a ground based sample of faint galaxies to / = 25 is peaked at rhi — 0.3 arcsec, 
in good agreement with that found at similar magnitudes in the HST medium deep survey. 

A morphological selection based on the asymmetry and concentration parameters gives a distribution 
similar to the one found by Abraham et al. (1996) from HST data. From this distribution a sample 
of normal spirals has been selected. Photometric redshifts have been obtained for this sample using the 
available multicolor catalog. 

A clear distinction in the morphological properties of galaxies fainter or bighter than Mb ~ — 19 has 
been found, especially in the 0.4 < z < 0.8 redshift interval. In particular the fainter galaxies show a marked 
concentration toward small (r^ < 3 kpc) sizes (see Fig. 5). 

The comparison of the observed galaxy sizes with those predicted by hierarchical clustering CDM models 
at intermediate redshifts (z ~ 0.5 — 1) indicates an excess of faint, small-size galaxies with respect to the 
model predictions (see Figs. 6, 8, 11). A better agreement is found at brighter (Mb < —19) magnitudes. 
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Different reasons could be responsible for this excess: 

a) The ad hoc assumption that all the disks have masses and angular momenta that are a fixed fraction 
of those of their DM counterparts (jd ~ rrid — const.) is incorrect for small galaxies. In fact, numerical 
N-body simulations indicate that the gas loses most of its angular momentum to the DM during galaxy 
assembly (Navarro & Benz 1991; Navarro & White 1994; Navarro & Steinmetz 1997), although the effect of 
Supernovae feedback (not included in the simulation) is still unclear. In this case, however, the problem of 
explaining why the assumption j^ « md ~ 0.3 works for more massive galaxies should be solved. 

b) The faint-end of the galaxy luminosity distribution evolves slower than predicted by hierarchical 
models. The observed population of small galaxies corresponds to the turn-off in the galaxy luminosity 
function at the faint end, and those galaxies constitute a population evolving in a physically different way 
from the body of the distribution due, for example, to the lack of aggregation of objects of small circular 
velocity. This would break the self-similarity of the hierarchical clustering evolution at small galactic scales. 

c) The luminosity of the satellite small galaxies within a given DM halo does not take into account star 
bursting phenomena due to merging interactions within the halo. For a fixed galaxy luminosity this implies 
an excess of small size objects. 

d) The evolution of the number of galaxies within a halo is poorly described by semianalytical models. 
The two body aggregation processes in fact could produce a strong number evolution resulting in an excess 
of small objects as the redshift increases. 

These results have to be checked on a larger sample where the morphological parameters are measured 
on high resolution images. The VLT telescopes could be ideal in this respect in terms of field of view, 
photometric depth and optical quality. Moreover, a caution should be emphasized when we associate the 
galaxy morphological properties derived from the observations to the unknown dynamical properties of 
the faint galaxies. All the previous arguments apply only if these faint objects are sustained by rotation. 
Rotational curves of galaxies with / ~ 25 can not be measured with the current technology on 8m class 
telescopes. Thus, the statistical approach presented in this paper represents a first step to include the galaxy 
morphology in the CDM models as a new constraint for galaxy formation and evolution. 

We thank the referee Dr. S. Odewahn for his helpful comments. We are grateful to S. Zaggia for 
providing HST archive images and A. Renzini for discussions. The ultraviolet observations of the NTT 
Deep Field were performed in SUSI-2 guaranteed time of the Observatory of Rome in the framework of the 
ESO-Rome Observatory agreement for this instrument. 



A. The Model of Galaxy Formation 

The model proceeds through the following basic steps: 

a) We define a time-circular velocity grid, with grid size At and Avc- Since a given DM halo may 
contain several galaxies, we denote with Vc the circular velocity of the halos (whose DM mass is M) and 
with Vc that of the galaxies; the latter represents the depth of the galactic potential wells corresponding to 
a mass m. 

The number density N{Vc,t{zi)) of Dark Matter (DM) halos with circular velocity Vc present at the 
initial redshift Zi is taken from the Press & Schechter expression. Initially, we associate to each DM halo 
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a galaxy whose circular velocity Vc is the same of the halo. To each Vc, a halo mass M = V^ / [10 G H{z)] 
and a virial temperature T^ are associated. The initial gas content associated to all the halos corresponds 
to the universal baryonic density Qf,. The initial stellar content is taken to be 0. For each circular velocity, 
the baryonic content is divided into a cold gas phase (the one able to cool within the halo survival time r/), 
with mass mc, into a star phase, with mass to*, and the complementary hot phase, with mass tti/j. 

b) For each galaxy circular velocity Vc, we compute the variations of the mass of the gas components 
and of the stellar content due to merging of DM halos, possibly followed by galaxy merging. At each time 
step, we first compute i) the probability P{V^, Vc, t) that a halo at time t with velocity FJ is included in 
a halo with velocity Vc at time t + At, and ii) the associated halo merging rate dP{Vc),t\Vc)/dVcdt: these 
are computed in the framework of the so-called Extended Press & Schechter Theory (EPST, see Bond ct al. 
1991; White & Frenk 1991; Bower 1991; Lacey & Cole 1993). 

Then we compute the number of galaxies Ng{vc, Vc, t) with circular velocity Vc in a halo with velocity 
Vc, implementing the canonical recipes of semi-analytical models. According to these, during halo merging, 
a galaxy contained in one of the parent halos contributes to enrich the dominant galaxy of the common 
halo (of velocity Vc) if its coalescence time Tm{vc) is less than the halo survival time ti{Vc) (predicted in the 
EPST). In the remaining case (t„i > rj) the galaxy retains its identity. 

Thus, for galaxies with Vc = Vc, the positive increment in the number of galaxies of given Vc in halos 
with Vc is linked to the number increment of DM halos with Vc [whose contribution from progenitors VJ is 
N{Vc, t) d^P{Vc, t\Vc)/dVcdi\ through the probability profe[Tm(w^) < T((Vc)] of forming one dominant galaxy 
merging those contained in the parent halos; for galaxies with Vc =/= Vc the positive increment is linked to the 
complementary probability l—prob\Tm{vc) < ti{Vc)\ ■ The overall negative increment is due to disappearence 
of satellite galaxies of velocity Vc in halos Vc , when the latter merge to forme larger halos V^ . Thus 

Ng{vc, Vc, t + At)^ Ngivc, Vc, t) 



Jv, dVcdt 

where 5{vc, Vc) = 1 for Vc — Vc and elsewhere. 

NgT{Vc) = / dv'c Ng{v'c, Vc, t) is the total number of galaxies in the halo VJ. The merging time of the 
galaxies is parametrized (in terms of the DM masses M{Vc) and M'(T/,')) as t„i = To{M' /M)"'"^'"'!' according 
to Cole et al. (1994), while the probability prob[Tm{v'c) < ti{Vc)] for a DM halo of velocty Vc to have a 
survival time t; larger than a given value has been computed by Lacey & Cole (1993) (see their eq. 2.21) 
in the framework of the EPST. The initial condition is such that in each halo there is only one galaxy (i.e., 
Ng{Vc,Vc,to) ~ N{Vc,to) for Vc — Vc and elsewhere). From the ratio Ng{vc,Vc,t)/N{Vc,t) we compute 
the probability f{vc, Vc, t) of finding a galaxy with circular velocities Vc in a DM halo of circular velocity Vc- 
Then the total number density of galaxies with circular velocity Vc (and the corresponding probability density 
p{vc, t)) is computed from Ng{vc, t) = J dVc N{Vc, t) f{vc, Vc, t), where the halo distribution N(Vc, t) takes 
the canonical Press & Schechter form. 

Finally, for each Vc we compute the probability distribution p{v'^,Vc,t) of having a galaxy at time 
t with circular velocity w^ included in a galaxy with circular velocity Vc at time t + At. This is given 
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by p{v'^,Vc,t) = J^ J^ dVcdV^ P{V^,Vc,t) f{v'^,V^,t) f{vc,Vc,t); from this we also derive the associated 
conditional probability p{v'^,t\vc) ~ p{v[.,t\vc)/p{vc,t) that, given a galaxy with velocity Vc at time i, a 
galaxy v'^ has been included in it at t — At. 

The average increments of the cold gas mass and the stellar mass in galaxies with circular velocity Vc 
are given by the equations in White & Frenk (1991); in differential form these read 

mc{vc,t + At) = mc{vc,t) + At I dv'^^fp^ p{v'^,t\vc)mc{v'^,t) (A2) 

m*(«c,t + At)=m,(wc,t) + At / dv'^-fp^ p{v'^,t\v^)m4v'^,t) . (A3) 

Jo Ng(Vc,t) 

The hot gas reservoirs are assumed to be retained in the DM halos rather than in individual galaxies. The 
amount of hot gas in halos with velocity Vc is given by 

m,(K,i + At) = m.(K,t) + At^''^dK'^|^^^^|^m.(T/',t). (A4) 

c) The mass of the cold phase is further incremented by the cooling: the gas that cools in a halo with Vc 
in the time step At is AmdVc) — 47rr^Q^j Pg{rcooi) ^fcooh where Pg is the gas density and the dependencies 
of the coohng radius rcooi on Vc and t are given, e.g., in Somerville & Primack (1998). Such cooled gas is 
assigned to the galaxies which have the same circular velocities of their DM halo, so that the average hot 
gas in galaxies with circular velocity Vc is Amdvc) = Amc{Vc) Ng{Vc, Vc)/ J dV/ Ng{Vc, V/)- 

The amount of cold gas which forms stars is then computed and the corresponding mass is transferred 
from the cold phase to the star phase. Such amount is Ato*(wc, t) — M* At and it is regulated by the star 
formation rate m*(wc,t) = [mc(uc,t) — m*(wc,t) — TO/i(wc,t)]/r*, which as in Cole et al. (1994). The timescale 
Tt{vc) = T° (wc/300 km/s)"' is parametrized according to Cole et al. (1994). 

A final transfer from the cold phase back to the hot phase, due to supernovae feedback, involves 
the amount Amh{vc,t) — jiiv^ Am^Vc^t) where /3 = (wc/^h)"'' is defined by the free parameters v^ 
and ah- The total hot gas returned to a halo with velocity Vc by the member galaxies is Avrihiyc) = 
/p" dvcNg{vc,Vc,t)Amh{vc,t). 

d) For each circular velocity Vc, the associated luminosity at wavelength A is computed by the convolution 

^A = (l/r) / Mt-t')mdt')dt' , (A5) 



^0 

where the spectral energy distribution of luminous stars (px (t) is taken from the stellar population synthesis 
models by Bruzual & Chariot (1996). The output normalization depends on the free parameter F, i.e., 
the fraction of "stellar" mass locked in non- luminous objects (planet-like objects). The dust absorption is 
included in the computation with the same procedure described in Somerville & Primack (1998). 

f ) We increment the current time by At and repeat the whole procedure from step b) until the output 
time is reached. 

As explained in the text, we choose our parameters in order to fit both the local observed Tully-Fisher 
relation and the local observed B-band luminosity function. For the star formation parameters, the above 
procedure yields the following set: 
T^ — 2 Gyr, a^ = —0.5, Vh = 100 km/s and q;/i=5 for the stellar formation. 
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r = 1 and Qh = 0.06 for the normalization of the niass-to-luniinosity ratio. 

To = 0.1 dynamical times, and «,„ = 0.2 for the galaxy merging rate. In addition, a cutoff in the cooling for 

halos with Vc > 500 km/s is included (see Somerville & Primack 1998). 
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FIGURE CAPTIONS 



Fig. 1 The photometric redshift distribution of galaxies with / < 25 in the NDF. 

Fig. 2 Deconvolved half light radii as a function of true values in simulated data. Filled points are for the 
bright objects, empty for the faint ones. Error bars are one sigma confidence intervals. 

Fig. 3 Half light radii obtained deconvolving the HST and NDF images of the same objects in the I band. 
The stright line represents equal values. 

Fig. 4 Distribution of the half light radius for galaxies with / < 25. The values are in log (arcsec). 

Fig. 5 Disk scale length (r^ = r/ij/1.68) as a function of luminosity for galaxies with / < 25 in two redshift 
intervals. Filled circles correspond to the disk-dominated subsample whose selection is shown in Fig. 6 and 
discussed in Sect. 4. 

Fig. 6 Histogram of the differential size distribution for two redshift intervals. The plotted numbers are 
galaxies per 5 arcmin^ field per kpc. 

Fig. 7 Distribution of the galaxies with / < 25 in the asymmetry - concentration plane (see text). The two 
straight lines delimit regions corresponding to different spectral types according to Abraham et al. (1996). 
The region inside the lines is occupied by normal spirals. The outside region towards low concentration and 
high asymmetry is occupied by galaxies with irregular morphology. Finally, the outside region towards high 
concentration and lower asymmetry is occupied by SO and few elliptical galaxies. 

Fig. 8 Size distribution of galaxies as shown in Fig. 6 but for galaxies having compactness and asymmetry 
typical of normal spirals. These galaxies are inside the region delimited by the two straight lines in Fig. 7. 
The continuous curve is the distribution predicted by the CDM model described in Sect. 5. To compare with 
the observed data in our sample, we have considered in the predicted counts a magnitude limit at / < 25. 
This limit has the effect of selecting brighter (and hence larger) objects at higher z. Such an effect tends to 
balance the decrease with z of the disk size predicted by CDM models (as shown in fig. 10 at given absolute 
magnitude); thus the predicted distributions in the two redshift intervals show an apparent lack of evolution 
with z. 

Fig. 9 Toj): TuUy-Fisher relation; the shaded area represents the region allowed by the observations; the 
continuous curve is the predicted T-F relation for the dark matter circular velocity while the dashed curve is 
the T-F relation for the disk circular velocity. Middle: the predicted local luminosity function; the data are 
adapted from Zucca et al. (1998). Bottom: the predicted total star formation rate as a function of redshift. 

Fig. 10 The fractional size distribution of galaxies predicted by the CDM model at different z for two ranges 
of absolute magnitude. The distribution results from the distribution of angular momentum A, which is 
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related to the disk radius Rd and circular velocity Vd by eq. 5. The dependence on the circular velocity 
is translated into a dependence on absolute magnitude through the TuUy-Fisher relation predicted by the 
model and shown in the first panel of Fig. 9. Note the growth of the average disk size with cosmic time 
typical of CDM scenarios. 

Fig. 11 Distribution of galaxies in the luminosity-size plane in four redshift intervals. The disk radii are in 
kpc. The shaded region represents the region allowed by the model. The upper and lower lines correspond 
to the 90% and 10% points of the angular momentum distribution. The solid line corresponds to the 50% 
point of the same distribution (see the text for details). 
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